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Abstract 

The  incomplete  beta  function  arises  in  various  statistical 
problems.  It  is  known,  for  example,  that  the  tail  probability  of  the 
binomial  distribution  can  be  expressed  as  an  incomplete  beta  function 
This  paper  gives  some  results  on  a  monotonicitv  property  of  the 
incomplete  beta  function.  The  given  results  are  shewn  to  have  appli¬ 
cation  in  a  problem  of  ranking  and  selection. 
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1.  Main  results.  There  are  given  below  two  theorems  on  the 
monotonicity  property  of  an  incomplete  beta  function.  An  applica¬ 
tion  of  the  given  result  is  shown  in  the  next  section.  Let 


Ip (a,b) 


r (a+b) 
r  (a)  r  (b) 


p  a-1  b-1 

x  (1-x)  dx, 

0 


a  >  0,  b  >  0,  0<p<l 

denote  the  incomplete  beta  function.  Let  a  =  £n+y  and  b  =  (l-£)n+$, 
where  0  <  £  >  1,  n>0,  y  j>  0 ,  5  :>  0.  The  following  theorem  estab¬ 
lishes  a  monotonicity  property  of  I  ($n+y,  (l-^)n+c)  in  terms  of 

P 

n.  Let  n'  >  n  and 

f (p)  =  I  (n'^+y,  n'  (I-*)  +5)  -  I  (n£+y ,  n(l-^)  +i). 

P  P 

Theorem  1.1.  Let  0<p  <1.  Iff(p)  <  (>)  0  then  f(o)  <  (>)0 

- ro  o  —  —  —  — 

for  p  <_  (■>)  pQ. 

Proof:  Let 

n^+y-1  n(l-*)+5-l 

gn (x)  =  x  (1-x)  /  B (n^+y ,  n(l-£)+c) 

denote  the  beta  density  function,  and  let 
h (x)  =  gn, (x)  /  gn  (x) . 

Clearly,  h(x)  is  nondecreasing  (nonincreasing)  in  x  for  x  <_  -"(_>)  i . 
Consider  the  function 

'P 

f(p)  =  !  (h(x)-l)g  (x)  dx. 

I  r\ 


(1.1) 


(2) 


As  x  varies  from  0  to  1 ,  the  integrand  on  the  right  hand  side  of 
(i.i)  is  either  negative  throughout  or  it  changes  sign  from  nega¬ 
tive  to  positive  and  then  to  negative.  Since  f ( 1 )  =  0,  it  follows 
that  f(p)  changes  sign  once  from  negative  to  positive  as  p  varies 
from  0  to  1.  The  conclusion  of  the  theorem  follows  immediately.  [] 


Let y =  0,3=1.  If  n  and  n^  are  integer  valued  then 


I  (ru,  n(l-U+l)  =  7  (^)pr(l-p)n~r 

p  r  =  ;n  r 


represents  the  tail  probability  of  a  binomial  distribution.  Let 
*<j(n)  =  a  log  7(n)/3n  denote  the  digamma  function,  and  let 

A  (n)  =  -1  (n)  -  ;>|n;)  -  (l-£ )  v  (n  (1--;) )  +  Uog?  +  (1-;)  log  ( 1-:) -1/ 

Using  the  integral  formula  for  the  digamma  function,  given  by 


1  f x  ?  t  t  -12  2-1 

■An)  -  log  n  -  - - 2  t(e*  -1)  (t  +n  )  dt 

‘a  •  0 


we  have 


(1.3)  Am)  =2  c(e2  1 )  ~  ^  p +  -y  L-3 - -  -  dt 

•a  t  +n“-;  t  +n“(l-l) "  t“+n- 


0  t2+n2-:2  t2+n2(l~:)2  t2+n2 


=  0  . 


Consider  the  function 


(3) 


C  (x)  =  ^ log  x  +  (l-$)  log  (1-x) --;iog^  -  {!-£) log ( 1-i)  +  i 


0  <  x  <  1. 

Clearly  C(x)  is  a  concave  function  of  x.  Let  pp  and  qp  denote  the 

roots  of  the  equation  C(x)  =  0,  where  p  <  %  <  q  .  Note  that 

n  n 

p  and  q  -  l.  as  n  <=.  We  have 
*n  ^n 


(1.4) 


32Ip(^n,  (l-On+D/apan  =  f-(~  { +Tf 


Uog  p 


(l-C)log(l-p)  +  i  (n+l)  -  CMnC)  -  ( 1-?)  v  (n  ( 1- ;)  +1)  ] 


n£-l.,  ,n(l-') 

_  P  (1-P) 

B(n5,n(l-- >+l) 


[C  (p)  +  A  (n)  ]  . 


The  second  equality  in  (1.4)  follows  from  the  relation  ;(n+l)  = 
i'(n)  +  1  In  view  of  (1.3)  we  have  that  the  right  hand  side  of 

(1.4)  is  negative  for  p  <  p  and  p  >  q  . 

3  *•  —  c  n  *■  —  *n 

Since  31  (n^  ,  n  ( 1- 1 )  +1) /3n  -*•  0  as  p  -*■  0  ,  it  follows  that 

ir' 

(1.5)  51  (n^,n(l-^)+15/5n<0 

P 

for  p  <  p  •  Using  the  relation  I  (n£ , n (1-^) +1)  =1-1.  (n(l-^)+l, 

"  —  ^  n  ^  p  1-p 

we  find  that  the  reverse  inequality  holds  in  v  1 . 5 )  for  p  q  q  .  We 
have  proved  the  following  result. 


Theorem  1.2.  The  incomolete  beta  function  I  (n£ , n ( 1- ; ) +1 )  is 

-  -  p 


(4) 


decreasing  in  n  for  p  <_  and  increasing  in  n  for  p  ->  q^,  where 

p  and  q  are  the  roots  of  the  equation  C(x)  -  0. 
n  n 


Remark  1.1.  It  can  be  shown  that  the  monotonicity  property 

given  by  Theorem  1.2, holds  also  for  the  function  I  (n£,n(l-£)), 

P 

where 

n  , 

I  (n;,n(l-;))»  l  (n"  )pr(l-p)n 

p  r  =  c,n 

if  n  and  n^  are  positive  integers. 


Remark  1.2.  From  Formula  (le.6.2)  of  Rao  (1966)  we  have 


C{X)  <  n  '  (“  2^~  ' 


Therefore,  0  < 


and  0  <  q 
^n 


2.  Application.  Consider  the  following  problem  of  ranking  and 

selection.  There  are  given  k  populations  with  cumulative  distribution 

functions  (cdf)  F^(x)  =  F^  (i=l,...,k),  and  a  number  a  with  0  <  a  <  1. 

The  distribution  functions  are  unknown  but  they  are  assumed  to  be 

continuous.  Let  ^  denote  the  a-quantile  of  F^ .  It  is  assumed  for 

simplicity  that  -1^  is  uniquely  determined  for  each  i  =  l,...,k. 

Given  a  sample  of  n  observations  from  each  population,  it  is  required 

to  select  the  population  associated  with  the  largest  value  of 

called  the  "best"  population.  We  snail  assume  that  nu  is  integer  valued. 

Let  x.  ■  denote  the  ith  order  statistic  in  the  sample  from  F .  ,  and 
l]  J  i 

let  j  =  nn.  Suppose  that  the  population  associated  with  the  largest 
value  of  x. ■  is  selected  as  the  best  oooulation. 

lj  ‘  * 


Let  the  ith  population 


(5) 


be  the  best  population.  Then  the  probability  of  a  correct  selection 
(PCS)  is  given  by 

(2.1)  PCS  =  n("“J)  uj"1(l-u)n"j  £=1 Ip  F-l(u) (na,n(l-a)+l)du. 

'0  t=j=i  t  i 

For  large  n,  the  right  hand  side  of  (2*1)  is  approximately  given  by 

k 

(0  n  1  ~  (na,n  (l-a) +1)  . 

t=l  F  (5?) 
tfi 


If  it  is  assumed  that  the  a-quantile  of  the  best  population  is 
sufficiently  larger  than  the  a-quantile  of  each  of  the  remaining 
populations,  in  the  sense  that 

Ft  »  t  f  i 

where  is  a  given  positive  number,  then  (2.2)  is  minimized  for 
F  ( - a)  =  a  +  £ ,  tfi. 

L.  J. 

Therefore,  the  minimum  probability  of  a  correct  selection  is  approx¬ 
imately  given  by 


(2.3)  min  PCS  =  (I  (na ,n ( i-a) +1) ) k_i 


2 

The  right  hand  side  of  (2.2)  is  increasing  in  n  tor  s  >  Z  by  Theorem  1-2 

V 

and  Remark  1.2.  Thus  a  minimum  value  of  n  can  be  determined  for 


the  given  selection  problem,  for  which  the  probability  of  a  correct 
selection  is  at  least  as  large  as  a  given  number  ?*  (4  P*  ^  1) . 


(6) 


The  problem  of  selecting  the  best  population  for  the  largest 
1-quantile, has  been  considered  by  Rizvi  andSobel  (1967)  and  Sobel 
(1967)  . 

In  the  application  given  above,  the  problem  of  selecting  the 
population  associated  with  largest  median  value,  that  is,  when 
a,  =  — ,  is  of  special  interest.  For  this  case  the  quantity  inside 
the  square  bracket  on  the  right  hand  side  of  (1.4)  is  given  by 


(2.4) 


C(p)  +  A(n) 


p (n)  -  y (|)  +  ~  log  (p(l-p) ) 
-(yt2—)  -  'p(n)  +  log  ( 4p  ( 1-p)  )  )  . 


Let  pQ  and  l-pQ  be  the  values  of  p  obtained  by  equating  the  right 
hand  side  of  (2.4)  to  zero.  Then  I  (~,  ~  +  1)  is  decreasing  in  r.  for 

p  2  2 

p  <  n  and  increasing  in  r.  for  p  >  1-d  .  To  illustrate  our  result, we 
r  -  o  o 

give  below  values  of  I  (§,  ^  +  1)  for  p  =  (.5,  .55)  and 

p  2  2 

n  =  2,  4,  8,  (2)14,  20,  100.  It  appears  from  the  table  that 

Il/2'§'  §  +  D  decreases  as  n  varies  from  2  to  10  and  increases 
thereafter . 


2 

4 

8 

10 

12 

14 

20 

.50 

.7  500 

.  6875 

.  6367 

.  6230 

.6123 

.6047 

.  5881 

.  55 

.  7975 

.  7585 

.  7396 

.  7384 

.  7393 

.7414 

.7505 
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